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We show that non-interacting disordered electrons on a Bethe lattice display a new intermediate 
phase which is delocalized but non-ergodic, i.e. it is characterized by Poisson instead of GOE 
statistics. The physical signature of this phase is a very heterogenous transport that proceeds over 
a few disorder dependent paths only. We show that the transition to the usual ergodic delocalized 
phase, which takes place for a disorder strength smaller than the one leading to the localization 
transition, is related to the freezing-glass transition of directed polymers in random media. The 
numerical study of level and eigenstate statistics, and of the singular properties of the probability 
distribution of the local density of states all support the existence of this new intermediate phase. 
Our results suggest that the localization transition may change nature in high dimensional systems. 



After more than a half century of Anderson localiza- 
tion, the subject is still very much alive [l[ as proved 
by the recent observations of Anderson localization of 
atomic gases in one dimension Q and of classical sound 
elastic waves in three dimensions Also on the the- 
oretical side several questions remain open: Although 
there is by now a good understanding of the localiza- 
tion transition in low dimensional systems, culminating 
in a functional renormalization group analysis by a 2 + e 
expansion jij, the behavior in high dimensions, in par- 
ticular the existence of an upper critical dimension and 
the relationship with Bethe lattice analysis [B[ , is still an 
issue. Recently, there has been a renewal of interest on 
this problem because of its relationship with Many-Body 
localization, a fascinating phenomenon that should take 
place for disordered isolated interacting quantum systems 
and lead to a new kind of phase transition between a 
low temperature non-ergodic phase — a purely quantum 
glass — and a high temperature ergodic phase (a metal 
in the case of disordered electrons) la]. A paradigmatic 
representation of this transition [g, |7| is indeed Ander- 
son localization on a very high dimensional lattice — the 
Fock space — which for spinless electrons consists in an N- 
dimensional hypercube (N ^> 1 is the number of sites of 
the lattice system). Localization had an impact on sev- 
eral fields, in particular Random Matrices and Quantum 
Chaos. As a matter of fact, in the delocalized phase the 
level statistics is described by random matrix theory and 
generally corresponds to the Gaussian Orthogonal En- 
semble (GOE), whereas instead in the localized phase is 
determined by Poisson statistics because wave-functions 
close in energy are exponentially localized on very dis- 
tant sites and hence do not overlap; thus, contrary to the 
GOE case, there is no level-repulsion and eigen-energies 
are distributed similarly to random points thrown on a 
line. The relationship with quantum chaos goes back to 



the Bohigas-Giannoni-Schmidt conjecture, which states 
that the level statistics of chaotic (or ergodic) systems is 
given by random matrix theory, whereas integrable sys- 
tems instead are characterized by Poisson statistics [8|. 
This result can be full y w orked out and understood in the 
semi-classical limit [9l.llu] : for a quantum chaotic system, 
in the Ti — > limit, wave-functions at a given energy be- 
come uniformly spread over the micro-canonical hyper- 
surface of the configuration space. They are fully delocal- 
ized as expected for an ergodic classical system that cov- 
ers regions with same energy uniformly. Instead, quan- 
tum non-ergodic models, such as integrable systems, are 
characterized by Poisson statistics and localized wave- 
functions. All those results support a general relationship 
between delocalization-GOE statistics-ergodicity (simi- 
larly between localization-Poisson statistics-lack of er- 
godicity). In this work, following up a suggestion in Q, 
we show that surprisingly non-interacting disordered 
electrons on a Bethe lattice instead display a new in- 
termediate phase, which is delocalized and yet still not 
ergodic (henceforth "not ergodic" will be a synonym of 
"characterized by Poisson statistics"). Since the Bethe 
lattice should be representative of high dimensional lat- 
tices, our results suggest that for large enough dimen- 
sions the localization transition may change nature and 
be characterized by a new intermediate phase, which ac- 
tually could also play an important role in Many-Body 
localization [llj . 

The model we focus on consists in non-interacting spin- 
less electrons in a disordered potential: 
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where the first sum runs over all the nearest neighbors 
couples of the Bethe lattice, the second sum runs over 
all N sites; cj, c; are fcrmionic creation and annihila- 



tion operators, and t is the hopping kinetic energy scale, 
which we take equal to 1. The on-site energies €4 are 
i.i.d. random variable uniformly distributed in the inter- 
val [-W/2, W/2}. The Bethe lattice is defined as a k+ 1- 
random-regular-graph, i.e. a graph chosen uniformly at 
random among all graphs of TV sites where each of the 
sites has connectivity k + 1 [13. This lattice is like a 
Cayley tree wrapped onto itself [13| . 
Localization on the Bethe lattice was first studied by 
Abou-Chacra, Anderson and Thouless [H and then later 



by many others, see [1J-|20( and refs. therein. Many sim- 
ilarities, but also few important differences, with the 3d 
behavior have been found. They mainly concern the criti- 
cal properties, which display exponential instead of power 
law singularities [3, EH, an d the Inverse Participation 
Ratio, defined as IPR = Yli=i IV^I 4 (V'i i s tne value of the 
wave- function on site i), which was conjectured to have a 
discontinuous jump at the transition from aO(l) toward 
a 1/N scaling [TBI] . The level statistics was conjectured 
to display a transition from GOE to Poisson concomi- 
tant with the localization transition. However, numeri- 
cal studies didn't fully support this claim 17[. Moreover, 
arguments put forward in [7J indicated that the two tran- 
sitions could actually not coincide. 

In the following we indeed unveil that Anderson local- 
ization on the Bethe lattice is more subtle than origi- 
nally thought and that some of the previous conjectures 
have actually to be revised. We focus on the k = 2 case 
(connectivity three) and on the middle of the spectrum, 
E = 0. Previous works determined by studying transmis- 
sion properties and dissipation propagation that the lo- 
calization transition takes place at Wl — 17.5 (Bl flil lioj] . 
In order to analyze the level statistics and clarify its rela- 
tionship with the localization transition, we have diago- 
nalized the Hamiltonian for several system sizes N = 2", 
with n = 9,10,11,12,13. For each N, we have aver- 
aged over both the on-site quenched disorder, and on 
random graph realizations (2 20_n samples). Since we 
are interested in E = 0, we only focused on 1/8 of the 
cigenstates centered around the middle of the band (tak- 
ing 1/16 or 1/32 of the states do not affect the results). 
We have studied the statistics of the level spacings of 
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En is the energy of the n-th eigenstate in the sample 
a. In order to avoid problems related to the unfold- 
ing of the spectrum, we measured the ratio of adjacent 
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gaps [2i|J: rn' = min{<s£ "'.^J/max^r,^}, and 
obtained the (sample averaged) probability distribution 
P(r), which is expected to display a universal form de- 
pending on the level statistics [2l[. We did find that 
P(r) converges to its GOE and Poisson counterparts for 
W < Wt and W > Wt respectively [Hf (the evolu- 
tion with the system size of the entire distribution P(r) 
is presented in the EPAPS). Surprisingly, however, Wt 
turns out to be definitely smaller than Wl- In order 
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FIG. 1. (r) as a function of W for different system sizes, 
showing that the transition from GOE to Poisson becomes 
sharper as N increases and all the curves cross at Wt — 14.5. 
Inset: Finite size scaling of the same curves showing data 
collapse ( P ~ 0.21). 

to determine precisely Wt we have focused on the finite 
size scaling behavior of (r). As shown in Fig. [TJ the de- 
pendence of (r) on W and N can be described in terms 
of a scaling function of the variable x = N^(W — Wt), 
where Wt — 14.5 and (3 ~ 0.21. For x — > —00 and 
x — > +00, we recover, as expected, the universal values 
( r )GOE — 0.53 and (r)p ~ 0.39 correspondingly to GOE 
and Poisson statistics. Fig. [T] clearly shows the existence 
of an intermediate phase for Wt < W < Wl- Inde- 
pendent evidence can be gathered by studying the be- 
havior of the average value of the IPR (averaged over 
samples and eigenstates close to E = 0). Only for 
W < 9 we do find the standard scaling as 1/N. Be- 
tween 9 < W < W T , we find that (IPR) - l/N a . The 
exponent a <G (0, 1] decreases (almost linearly) from 1 to 
and seems to vanish right at Wt- In the new interme- 
diate phase, Wt < W < Wl, the IPR is better described 
by a logarithmic fit, (IPR) ~ l/(log7V) b . The exponent 
b decreases and apparently vanishes for W f Wl- Hence, 
wave-functions appear to be delocalized in the interme- 
diate phase but on a number of sites so restricted that 
cigenstates are effectively independent, they do not re- 
pel each other and Poisson statistics holds. Complemen- 
tary results, supporting our findings, have been obtained 
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The previous analysis was performed for finite lattices. 
In order to overcome this limitation and obtain a nu- 
merically exact solution directly for an infinite lattice, 
we follow Abou-Chacra et al. [H] and analyze the recur- 
sive equations satisfied by the resolvent matrix G{z) = 
(H — zl)^ 1 on a Bethe lattice: 

G£- (E + ir,) = ei -E-ir,-J2 Gf-yi(E + ir,), (2) 

j'edi/j 

where di/j denotes the set of neighbors of i but j, and 
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FIG. 2. Power-law tails of the probability distributions of 
the imaginary part of Gi->j in the delocalized ergodic phase 
[W = 13, top panel), delocalized non-ergodic intermediate 
phase (W = 16, middle panel), and insulating localized phase 
(W = 18, bottom panel), showing the value of the exponent ji. 
Note that only the curves in the top panel show a convergences 
to a 77-independent value. 

is the diagonal element of the resolvent matrix 
on site i for a modified lattice where the edge between 
site i and j has been removed. A similar equation, 
where the sum on the LHS is over all neighbors of i 
allows one to get the local resolvent Gu and, from it, 
spectral properties such as the global density of states, 
p{E) = (l/N)J2 a 6(E-E a ) = lim^ (l/iV7r) £ 4 ImG lt , 
and (IPR) = lim^o rj\Gu{E + irj)\ 2 . Since the Gj_^-s are 
random variables, Eq. ([2]) leads to a self-consistent equa- 
tion on their probability distribution, P(Gi^j), which 
can be solved numerically using a population dynamics 
algorithm 0, [Til [20I ] . In our analysis we used a popu- 
lation of M = 2 26 elements. Note that the value of M 
needed to correctly sample the tails of the distributions 
increases when rj decreases, for M. = 2 26 we could only 
access rj in the range 10~ 2 — 10~ 6 . 

We found that in the localized phase, W > Wl, the 
probability distribution of the imaginary part of Gj_^j is 
singular in the limit 77 — > 0. Almost all values of ImG^j 
are of order rj except very rare ones. More precisely, in 
the small rj limit, P(ImGi^>j) has a scaling form f{x/rj)/rj 
for x ~ rj, where J dyf(y) = 1, and very fat tails: 
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FIG. 3. Exponent /i of the power law tails of the probability 
distribution of the imaginary part, of the real part, and of 
the modulus of Gi->j as a function of W, showing the two 
transitions at Wt and Wl- 

W < Wt we find that as rj — >• the whole P(G,_ ) . J -) has a 
non-singular ^-independent limit characterized by power 
law tails (with logarithmic corrections) with an exponent 
/i > 1, see the top panel of Fig. [5] (and the EPAPS for 
more details). All these findings are in agreement with 
previous results. What is surprising compared to previ- 
ous expectations is that there is not a direct transition be- 
tween these two regimes but instead, for Wt <W< Wl, 
we find an intermediate phase characterized by a mixed 
behavior of P(ImGi->j) and P(ReGi^j), that have a non- 
singular rj- independent limit [lil [l8| but, at the same 
time, display singular power law tails such as ((3]) (likely 
with additional logarithmic corrections) with an expo- 
nent stuck to 1 for the real part and changing from 1/2 
to 1 for the imaginary part, see the middle panel of 
Fig. [2] [24 |. Note that all the results we presented for 
G,:_>j also hold for Gu. The evolution of the exponent /x, 
which is reported in Fig. [3l clearly shows that a transi- 
tion, different from Anderson localization, takes place at 
W T . Since lim^ ImG, t = n S(E - £& a) )|V$| a one 
expects that the rare and very large values of ImG^ are 
due to similar rare events in \ip„J | 2 - This is confirmed 
by exact diagonalization of finite systems (see EPAPS): 
the distribution of \^^}\ 2 displays power law tails with 
the same exponent ji appearing in More results, 

including multi- fractal behavior, will be presented else- 
where [25j | . 

The nature of the transition taking place at Wt can be 
understood by a mapping to directed polymers in random 
media (DPRM) [2q |. In order to do this, we telescope the 
recursive equation on ImGi-^ and end up with the ex- 
pression [l8( | : 



with c being a constant O(l), and the exponent ji = 1/2 
in the whole localized phase. The power law behavior is 
cut-off for ImGi_>j ~ I/77, see bottom panel of Fig. [2] 
Note that since ji < 1 the tails give the leading contri- 
bution to the DOS, whereas the bulk part only a vanish- 
ing one. The probability distribution of the real part of 
Gi^j instead converges for rj —> to a stationary distri- 
bution with power law tails with exponent ji — 1. For where V are all the directed paths of length R of the 
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(rooted) Bcthe lattice originating from site i, and (i 1 , i") 
are all the edges, including i — >• j, belonging to the 
path V. (We neglected the ^-contribution since for 
Wt < W < Wl the typical values of ImG,^ are of 
the order of one whereas rj — > 0.) The sum in Q is 
over an exponential number of paths, k R ; in the large 
R limit, two cases are possible: the sum is dominated 
either by few paths that give a 0(1) contribution or by 
an exponential number of paths, k' R , each of them giv- 
ing a very small contribution but such that their sum is 
0(1) (k' is less than k and VF-dependcnt). Wc found 
that the transition between these two regimes happens 
precisely at Wt and it is related to the glass transi- 
tion of DPRM. Indeed, by introducing the edge-energy, 
Ui'-^i", by the equation e~ uli '^ i " = |GV->i" | 2 , one can re- 



interpret Y\ 



(i',i")£P 



as the partition func- 



tion for a DP on the Bethe lattice in presence of quenched 
bond disorder. This problem can be solved, even in the 
case of disorder correlated as the |GV-u" | 2 s [3, 0,0], 
by computing the generalized "free-energy" (also intro- 
duced in [ll|): 



d>(s) = lim lim — lo 

ri->0 iJ->oo Rs 




(5) 



It is easy to check that </>(s) is a convex function and that 
6(s = 2) = for W < W L since typical values of ImGj—^ 
are of the order of one. The average free energy of the 
directed polymer has a freezing-glass transition, called 
one-step Replica Symmetry Breaking, akin to the one of 
the Random Energy Model [26 



27] : by decreasing the 



"temperature" 1/s the generalized free energy decreases 
until the critical point s*, defined by d(j)(s)/ds\ s = s ± = 0, 
is reached; for s > s* the DP freezes and its free energy 
remains constant: (f>(s) = <fi(s*). In this glass phase the 
number of paths contributing to (|5|) is not exponential 
in R, as it happens for s < s*, but instead 0(1). It has 
been ri goro usly proved in [l8| (and indirectly found in Q , 
see [3, 25 1 for a detailed explanation) that the localiza- 
tion transition takes place for s* = 1, 4>(s*) = [28|. By 
computing <fi(s) using recursive equations we have indeed 
found that s* = 1 at W = Wl ; when diminishing W be- 
low Wl the value of s* increases and eventually reaches 
2 for W = W T - (see the EPAPS for details and plots of 
4>(s)). In consequence, the mapping on DPRM unveils 
that the intermediate phase corresponds to a regime in 
which there is derealization (s* > 1), but the sum in ^ 
is dominated by few paths only (s* < 2). The statisti- 
cal behavior of ImGi^j can also be understood in terms 
of DPRM: it was shown in [2(| that in the freezing-glass 
phase the DP partition function, and hence ImGj_>.j , dis- 
play power law tails with an exponent s*/2 and logarith- 
mic corrections; this explains the emergence of the ex- 
ponent \i that we indeed find equal to s*/2. Note, in 
particular, that s* = 1 and s* = 2 correctly match the 
values fi = 1/2 and fi = 1 at W = Wl and W = Wt- 



In conclusion we have shown that disordered non- 
interacting electrons on a Bcthe lattice display a new 
phase that is delocalized and non-ergodic. This comes to 
a large extent as a surprise for a model that was consid- 
ered to be exactly solved 0, Qj , 15 1 . Our results are not 
actually in contradiction with previous ones since those 
mostly focused on the localization transition at Wl ■ The 
only disagreement regards the behavior of the (IPR) that 
was su ppo sed to behave as 1/N in the whole delocalized 
phase |29|. This is likely due to a hypothesis of reg- 
ularity of the order parameter function intervening in 
the solution [l5[ that does not hold [3p| . The mapping 
to DPRM unveils the peculiar nature of the delocalized 
non-ergodic phase: the level broadening for a given site 
i, ImGn, is finite for 77 — >■ 0, but cutting an arbitrarly 
far away bond chosen on one of the few paths contribut- 
ing to dU substantially changes ImG^. Physically, this 
means that dissipation and transport are very heteroge- 
nous: particles can travel far away but on a few and very 
specific, disordered dependent, paths. Bethe lattice re- 
sults are often considered representative of the infinite di- 
mensional limit. In consequence our results suggest that 
for large enough dimensions the localization transition 
may change nature. This possibility, the implications 
for Many Body localization, the relationship with simi- 
lar intermediate phases [3ll433j and the connection with 
DPRM [Ull^] certainly are worth further investigations. 
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and results related to several points discussed in the main 
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FIG. 4. Probability distribution of the average gap ratio for 
different system sizes (from N = 512 to N = 4096) and for 
four different values of the disorder. The Poisson and GOE 
counterparts of P(r) are also shown. Top-left panel: Weak 
disorder, W = 8. The entire probability distribution P(r) is 
described by the GOE ensemble. Top-right panel: Moder- 
ately weak disorder, W = 12. P(r) is approaching the GOE 
distribution as the system size is increased. Bottom-left panel: 
Moderately strong disorder, W — 16. P(r) is converging to- 
wards the Poisson distribution as N is increased. Bottom- 
right panel: Strong disorder, W = 23. The entire probability 
distribution is described by the Poisson one, except at very 
low r, where convergence is exponentially slow due to finite 
size effects. 



Level statistics: 



The probability distribution of the 
gap ratio. 



In the case of Poisson statistics the probability distri- 
bution of the gap ratio r is given by P(r) = 2/(1 + r) 2 
(brown curve of Fig. 3|) , and its mean value is (r) p = 
2 log 2 — 1 ~ 0.386. For very strong disorder, in the lo- 
calized regime, the numerical results show that P(r) has 
indeed this form (except for very small values of r, for 
which the data approach the Poisson distribution expo- 
nentially slow due to finite size effects), as shown in the 
bottom-right panel of Fig. @]for W = 23. For interme- 
diate strong disorder, full convergence is not reached for 
our biggest system size. Nevertheless, it is clear that 
P(r) is converging toward the Poisson distribution as N 
is increased, as shown in the bottom-left of Fig. @] for 
W = 16. 

No analytical expression for the probability distribu- 
tion of the gap ratio is known in the case of GOE statis- 
tics but it can be easily determined numerically from ex- 
act diagonalization of large GOE random matrices (cyan 
curve of Fig. [4]). The mean value (r)aoE is equal to 
0.5295 ± 0.0006. Level repulsion in the GOE spectra 
manifests itself in the vanishing of the probability dis- 
tribution at r = 0. As expected, for weak disorder, the 
entire probability distribution P(r) is described by GOE, 
as shown in the top- left panel of Fig. 0] for W = 8. For 
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FIG. 5. W dependence of the exponents a and b describing 
the behavior of (IPR) as a function of the system size. 



intermediate low disorder, although full convergence is 
not reached for the biggest sistem size, P(r) is converg- 
ing toward the GOE distribution as N is increased, as 
shown in the top-right panel of Fig. 0] for W — 12. Right 
at Wt P(t) is described by a stationary iV- independent 
distribution (not shown), at least for the range of sys- 
tem sizes that we considered, which does not correspond 
neither to GOE nor to Poisson statistics and should be 
universal. 

These findings provide numerical evidence for the tran- 
sition in the level statistics from the GOE to the Poisson 
universality class. Such transition does not take place 
at the Anderson localization transition, Wl, as argued 
in [13], but at a smaller value of the disorder, Wt- Con- 
sider, for instance, what happens for W = 16 < Wl'- 
although the system is in the delocalized phase, P(r) is 
converging towards the Poisson distribution and (r) is 
approaching (r)p as N is increased (as we show by finite 
size scaling in the main text). 



Inverse Participation Ratio. 

The IPR of the n-th eigen-function of sample a is de- 
fined as IPR^ = J2iLi l^ni I 4 ' where ip^J is the value 
of eigen-state n of sample a on site i. In Fig. [5] we plot 
the values of the exponents a (red circles) and b (blue 
squares) describing the system size dependence of (IPR) 
(averaged over eigen-functions around E = and over 
all samples) as (IPR) - l/N a and (IPR) - l/(logiV) h 
respectively, as a function of the disorder. In the delo- 
calized ergodic phase, W < Wt, (IPR) ~ c/N a , with 
the exponent a decreasing from 1 at W ~ 9 to at 
W = Wt as the disorder is increased. In the interme- 
diate delocalized non-ergodic phase, Wt < W < Wl, 
the average (IPR) is better described by a logarithmic 
behavior, (IPR) ~ c/(log N) b . The exponent b decreases 




FIG. 6. Main panel: Q = prob(IPRl° !) > c/logJV) as a 
function of the disorder W for different system sizes, and for 
c = 1.5. Q grows from to 1 as the disorder is increased. 
The transition becomes sharper for bigger sizes. The curves 
obtained for different system sizes all cross around W ~ Wt- 
Inset: Finite size scaling plot of Q as a function of the variable 
AP (W — Wt) for different system sizes showing data collapse 
for Wt — 14.5 and 7 = 0.14 (the results should be indepen- 
dent of c for N large enough, however for the limited range 
of N at our disposal c ~ 1.5 provides the best crossing). 



as the disorder is increased and vanishes at Wl ■ Finally, 
for W > W L , (IPR) goes to a constant of 0(1) as TV is 
increased. 

In order to rationalize the unexpected behavior of the 
average IPR in the intermediate delocalized non-ergodic 

> 



(n) 



phase, we have analyzed the quantity Q = prob(IPR, 
c/logiV), defined as the probability that, for a given 
eigen-function n of a given sample a, IPR^*-* is bigger 
than c/ log N (c being a constant of order 1) averaged 
over eigen-states close to E = and over all samples. If 
the average IPR behaves as a power law as a function of 
the system size, (IPR) ~ 1/N a (as in the fully delocal- 
ized phase) then Q must vanish in the thermodynamic 
limit. On the contrary, if either (IPR) goes to a constant 
of O(l) for large sizes (as in the fully localized phase), 
or (IPR) ~ l/(logA^) 6 with b < l(as in the intermediate 
phase), then Q must approach 1 as TV — > 00. We indeed 
observe that Q grows from to 1 as the disorder is in- 
creased and that this transition becomes sharper as the 
system size grows. For a suitable choice of the constant 
c = 1.5, the curves obtained for different system sizes 
all cross at W ~ Wt and can be described in terms of 
a scaling function of the variable N 1 {W — Wt), where 
W T ^ 14.5 and 7 = 0.14, see Fig. [6l The value of 7 
is a bit different from the exponent used to collapse the 
(r) data. However, because of the rather limited range 
of system sizes and the small value of the exponents we 
do not consider this difference meaningful. Much larger 
system sizes would be necessary to get a precise estimate 
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FIG. 7. Full probability distributions of logImGi-_>j for dif- 
ferent values of the imaginary cut-off r) from 10" 2 to 10 -6 for 
four different values of the disorder. Top-left panel: W = 9. 
P(logTmGi->j) converges to a stationary 77-independent dis- 
tribution. Top- right panel: W = 13. P(loglmGi^j) is ap- 
proaching a stationary distributions although we are not able 
to observe full convergence for the smallest allowed values of 
77. The tails of the distributions are described by a power 
law (with logarithmic corrections) with an exponent [i = 1.3 
(black dashed line). Bottom-left panel: W = 16. In the 
intermediate delocalized non-ergodic phase P(logImG is 
bound to have a non-singular behavior for 77 — >• 0, as proven 
in [l^]. However, we are apparently still very far from con- 
vergence. The tails of the distributions are described by a 
singular power law (with logarithmic corrections) with an ex- 
ponent n — 0.76 (black dashed lines), whose coefficient be- 
haves as ?7 1_M . Bottom-right panel: W = 18. In the localized 
phase the limit 77 — > is singular. Almost all ImGi_>j are of 
order rj, except very rare ones described by a power law (with 
logarithmic corrections) with an exponent (j, = 1/2 (black 
dashed lines), whose coefficient vanishes as ^Jr\. 

of the exponents. 

Singular behavior of the local density of states. 

In Fig. [7] we show the full probability distributions of 
log(ImGj_>j) for several values of the cut-off 77. In the 
fully localized phase, W > Wl, the probability distri- 
bution of the imaginary part of Gi^j has a singular be- 
havior for 77 — > 0, as shown in the bottom-right panel of 
Fig. [7] for W = 18. Almost all the values of ImG^j are 
of order 77, except for very rare events (whose fraction 
vanishes as 77) , described by very fat power law tails with 
an exponent 1 + fj, = 1.5 (see Eq. ([3]) and Fig. [2]). Such 
tails give a contribution of 0(1) to the local density of 
states, whereas the bulk part only yields a vanishing con- 
tribution. The probability distribution of the real part 
of G;_>7 (not shown) converge to a stationary distribu- 
tion with power law tails with a VF-independent exponent 
1 + fj, = 2. 

For W < Wl, since the system is delocalized and 



the spectrum is absolutely continuous, P(lmGi^j) must 
have a non-singular limit as 77 — > 0, as shown rigor- 
ously in [18| . Indeed, for low enough disorder the whole 
P(log(ImGi_j.j)) converges for ?/ — > to a stationary non- 
singular 77-independent distribution, as shown in the top- 
left panel of Fig. [7] for W — 9. This is also the case for 
W = 13 (top-right panel of Fig. [7]), although we are not 
able to observe full convergence for the smallest value of 
77 we are allowed to consider (77 3> l/M). The tails of 
P(ImG;_i.7) are described by a power law (with logarith- 
mic corrections) with an exponent 1 + /i with fj, > 1. The 
probability distribution of the real part (not shown) also 
converges to a stationary distribution described by power 
law tails with the same exponent. 

Finally, in the bottom-left panel of Fig. [7] we show 
the behavior of P(log(ImGi^.,)) in the intermediate de- 
localized non-ergodic phase. The rigorous result of [l8[ 
guarantees us that the limit 7/ — > is not-singular and 
that a stationary distribution exists. However, for the 
smallest value of 77 that we are allowed to consider we do 
not observe any convergence yet. The most likely rea- 
son is that since typical values of ImGi-^ are finite (for 
?7 — > 0) but extremely small close to Wl, (they vanishes 
as exp(c/VM / L — W) [HI), values of 77 much smaller than 
the ones we can use are needed to observe convergence. 
We are instead able to observe in full detail that the 
distribution of ImGi^j displays singular power law tails 
(likely with logarithmic corrections) described by Eq. J3]) 
with an exponent 1/2 < /J, < 1, yielding a finite contri- 
bution to the local DOS. The probability distribution of 
the real part (not shown) converges to a stationary dis- 
tribution, whose tails are described by a power law (and 
logarithmic corrections) with a VF-independent exponent 
1 + n = 2 as in the fully localized phase. 



Wave-functions statistics in the delocalized phase 

In this section we report some results on the behav- 
ior of the wave-functions coefficients, l^ii | 2 > obtained 
from exact diagonalization of finite systems for W < Wl ■ 
Since lim^^ ImG,j = n J2 n ^( E ~ E n)\^i\ 2 it is natu- 
ral to expect that the rare and very large values of ImG;i 
are related to values of |^^| 2 much larger than the typ- 
ical ones (this is true only for W < Wl where the spec- 
trum is absolutely continuous, not for W > Wl where 
instead large values of ImGjj are produced by localized 
rare resonances). More precisely, we expect that power 
law tails for IitiGm should be related to power law tails 
for P(|^| 2 ). For strong enough disorder, we indeed 

observe that the tails of P(\ipn}\ 2 ) are described by a 
power law. The exponent value changes (increases) with 
the system size and approaches the value \x describing the 
tails of P(ImGii), as shown in Fig. [S] We interpret this 
finite size effect in terms of fluctuations of the effective 
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FIG. 8. Probability distribution of the wave functions am- 
plitudes, P(\ipni | 2 )' ^ or — 16 and different system sizes 
N = 2 n , with ?i = 9, . . . , 14, obtained from exact diagonal- 
ization of 2 18 ~" samples. The tails of the probability dis- 
tributions are described by a power law with an exponent 
1 + /x = 1.65 (black dashed lines) very close to the the value 
of the exponent of the tails of P(lmGi^j). The tails of 
PWnil 2 ) collapse on a single scaling function when plotted 
as a function of the scaling variable l^^i^N 1 ^ . 

value of Wl and the existence, for finite system sizes, of 
localized eigcnstate even for W < 17.5. These lead to 
spurious larger values of \if>„} | 2 that make the exponent 
of the tails effectively larger than its large-N limit. 

As previously shown in [l6[ we have found that \ip^nl 1 2 
display a multi-fractal behavior (more details will be pre- 
sented in a forthcoming paper (25[). Multi-fractality 
implies that the probability distributions of the wave- 
functions coefficients for different system sizes cannot be 
fully rescaled in terms of the variable \tp n " | 2 iV K using a 
unique exponent k. On the contrary, the value of k should 
depend on the range of \ip„J | 2 one focuses on. In partic- 
ular, the value K t describing the rescaling of the typical 
value of IV'i^l 2 is quite different from the ones, n es and 
K e i, describing the rescaling of extreme (very small and 
very big) values of the wave-function amplitudes. We 
have found both analytically and numerically (more de- 
tails will be presented in a forthcoming paper [25|) that 
in the intermediate phase the extreme tail of the distri- 
bution P(\ip n ,i\ 2 ) are given by 

const 



tf|iM 2(1+M) 

which means that K e i — l//i as shown in Fig. [8] for 
W = 16. There are logarithmic corrections to the previ- 
ous formula important to take into account to find that 

(IPR) = iV/odl^^PdVSl 2 )!^! 4 vanishes loga- 
rithmically in the intermediate phase. In agreement with 
the discussion above, we find that the whole curve cannot 
be rescaled using |V4"il 2 -^eZ- 



Computation of the "free energy" <J>(s). 

Before explaining the numerical computation of 4>(s 
we shall briefly discuss its physical meaning, see lg, [25 
[26l | for more details (in particular concerning the differ- 
ence between quenched and annealed averages). Denot- 
ing e~ R * s the contribution of a given path of length R, 
one can rewrite the sum ([5]) as an integral over all paths 
giving a contribution characterized by a value of / be- 
tween / and / + df times the number of such paths. By 
denoting the latter exp(i?S(/)), one ends up with the 
expression 



e n & 

V (i',i")eV 



dfexp(R[-sf + X(f)]) 



The value of / that dominates the integral for R — > oo 
depends on s. For small enough s, one finds that the 
saddle point value of f(s) is such that S(/(s)) > 0. In 
this regime an exponential number of paths, k' R (with 
k' = e s (*( s "), contributes to the sum. By increasing 
s, S(/(s)) decreases until the value s* is reached. At 
this point the generalized entropy E(/(s)) vanishes; in 
conclusion <fi(s) is related to the Legendre Transform of 
£(/), more precisely 4>'{s) = — S(/(s))/s 2 , and allows 
one to find out whether a finite or an exponential number 
of paths contributes to the sum in (0| . 

In order to compute <j>(s), on each edge of the lattice 
(and for each finite value of the imaginary regulator 77) 
we introduce the variable 



e n \°* 

V (i',i")eV 



(6) 



where V are all directed paths of length R originating 
from site i, and are all the directed edges (includ- 

ing i — > j) belonging to the path. It is straightforward to 
derive the following exact recursion relation for y^j : 



= \G 



E vr 

j'edi/j 



(7) 



where |Gj_j..,| s can be computed using Eq. (J5J). 
Eqs. ^ and (J7J) can be interpreted as an exact self- 
consistent equation for the joint probability distributions 
Q(Gi^j, yi^j), which can be computed by iteration using 
a population dynamics algorithm with arbitrary numer- 
ical precision. At each new i?-th step of the iteration, 
we can compute cf>(s, R, 77) as the average value of the 
logarithm of yt-yj over the distribution Q(G,_>j, Vi^.j), 
divided by Rs, as in Eq. ([7]). In order to take the 
R — > 00 limit, after about 10 4 iterations we extrapo- 
late the asymptotic value of <f>(s, R, 77) with a power law 
fit of the form: 4>(s,R,rj) ~ <j>{s,rj) + c(s,n)/R 1 ' (£ = 1 
for all values of W, s, and 77) where c(s, 77) is a con- 
stant which depends on s and 77. Finally, in order to 
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FIG. 9. Main panel: Free energy, 4>(s), as a function of s in 
the intermediate delocalized non-ergodic phase, for W = 15 
(blue circles), W = 16 (red squares), and W — 17 (green 
diamonds), showing the value of s* where d4>(s)/ds\ 3 = 3 * = 
and RSB occurs. Inset: Schematic sketch of the form of the 
wave- functions in the delocalized non-ergodic phase. Only 
few paths (the one colored in red) contribute significatively 
to the value of ImGi^j on the central site. 

take the limit 77 — >• 0, we computed <f>(s, rf) for several 
values of rj from 10~ 2 to 10 ~ 6 and used a fit of the form 
(j)(s, rf) = 4>(s) + Arf (we find that 6 ~ 1/4 for all values 
of W and s). Note that the limits 77 — >- and R — >• 00 
only commute when <f)(s) is finite: whenever <fi(s) = we 
find that the constant c(s,rj) diverges as r] — s- 0. 

At weak enough disorder, in the delocalized ergodic 
phase, W < Wt, <p(s) is a smooth decreasing function 
of s. It is positive for s < 2, decreases as s is increased 
and vanishes, with a positive derivative, at s = 2. Hence, 
in this case the contribution to ImGj-^, Eq. Q, comes 
from an exponential number of paths, since £(/(2)) > 
0. For W = Wt the glass transition takes place: the 
derivative of <j>(s) vanishes, i.e. s* = 2 (corresponding to 
fj, = 1, cfr. Fig. d and S(/(2)) = 0. In the intermediate 
delocalized non-ergodic phase, Wt < W < Wl, s* is 
found in the interval (1, 2): (f>(s) is positive for s < s*, it 
vanishes at s*, and remains zero for higher values of s. 
Correspondingly, the sum in |4]) is dominated by a finite 
number of paths only. In Fig. Owe show plots of 4>(s) in 
the intermediate delocalized non-ergodic phase, showing 
the value of s*, which coincides (within the numerical 
accuracy) with 2/j, (cfr. Fig. [3]). Note that in agreement 
with [l||, we did find s* = 1 at W L . 



Energy transport. 

It is very important to figure out what happens to 
energy transport and diffusion in the new intermediate 
delocalized-ergodic phase. In order to do that one can 
study the propagation of dissipation, i.e. of ImGj:_>j, on 
an open Bethe lattice where a thermal and particle reser- 
voir is attached to the boundary sites. An approximate 
way to perform this analysis corresponds to study the re- 
cursion equations, Eq. ([2]), where the imaginary regulator 
r] is set to a finite value of 0(1) only at the first step of 
the iteration, and rj = for all the successive iterations. 
Each iteration step corresponds to moving a generation 
further towards the bulk of the graph [B| . 

In the delocalized ergodic phase, W < Wt, after a cer- 
tain number of iterations P(lmGi^j) reaches a station- 
ary distribution which coincides with the //-independent 
limit shown in Figs. [5] and [71 as expected, and (ImG^ ) 
has a well-defined value in the bulk. 

In the fully localized phase, W > Wl, P(ImG^j) 
does not have a stationary distribution. As we move 
deeper towards the center of the graph, the whole distri- 
bution globally shifts exponentially towards smaller val- 
ues of ImGi-^j. Both (lxn.Gi-+j) and exp ((log ImG; _>.,)) 
go exponentially to as the number of iteration is in- 
creased. As expected, the system, even if coupled at the 
boundary to a reservoir, does not show any dissipation 
and decoherence in the bulk. 

In the delocalized non-ergodic phase, Wt <W< Wl , 
the situation is intermediate between the two cases de- 
scribed above. As exlained in the main text, the fact that 
the sum in (j4]) is dominated by few paths only implies 
that particle and energy propagation proceeds over few 
specific disorder dependent paths. Therefore, there is dis- 
sipation and decoherence in the bulk but it is extremely 
heterogeneous. Indeed, we find that the energy of the 
reservoir at the boundary penetrates towards the center 
of the system, since on almost all the sites of the bulk 
ImGi^j is finite, albeit possibly very small, as demon- 
strated by the fact that exp((logImG^j)) approaches a 
(very small) finite value. However, there are very rare 
resonances where ImGi-j.j is huge (possibly correspond- 
ing to the singular power law tails of P(ImG^j) found 
when rj is different from zero on all sites) . The behavior of 
such extreme events should be better characterized and 
understood. Qualitatively, we observe that the deeper we 
move towards the bulk of the lattice the rarer and bigger 
the resonances become. 



